Implementation of Grover's Quantum Search Algorithm 
in a Scalable System 
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We report the implementation of Grover's quantum search algorithm in the scalable system of 
trapped atomic ion quantum bits. Any one of four possible states of a two-qubit memory is marked, 
and following a single query of the search space, the marked element is successfully recovered with 
an average probability of 60(2)%. This exceeds the performance of any possible classical search 
algorithm, which can only succeed with a maximum average probability of 50%. 

PACS numbers: 



Quantum computers promise dramatic speedup over 
conventional computers in some applications owing to 
the power of entangled superpositions Q. Among the 
best-known quantum applications is Grover's search al- 
gorithm, which can search an unsorted database quadrat- 
ically faster than any known classical search [2j • A com- 
mon analogy for this searching algorithm is the prob- 
lem of finding a person's name in a phone book given 
only their phone number |3( : for N entries in the phone- 
book, this requires of order N queries. However, if the 
correlation between name and phone number is encoded 
with quantum bits, the name can be found after only 
about \/~N queries. While Grover's algorithm does not 
attain the exponential speedup of Shor's quantum factor- 
ing algorithm it may be more versatile, by providing 
quadratic gains for almost any quantum algorithm or 
accelerating NP-completc problems through exhaustive 
searches over possible solutions |y] . 

We implement the Grover search algorithm over a 
space of N=4 elements using two trapped atomic ion 
qubits Grover's algorithm has been implemented 

with ensembles of molecules using nuclear magnetic res- 
onance 0, 0, lllij with states of light using linear opti- 
cal techniques [l2Ml3j, and with Rydberg states within 
individual atoms |14j . None of these systems are scal- 
able however, as they require exponential resources as 
the number of qubits grows. The implementation of 
Grover's algorithm reported here complements the reper- 
toire of multi-qubit quantum algorithms recently demon- 
strated in the scalable system of trapped atomic ions 
(lil Uti Iril Hif . Unlike these earlier ion trap demon- 
strations, we use magnetically-insensitive "clock state" 
qubits and particular entangling gates that are uniquely 
suited to such qubits while remainin g in sensitive to ex- 
ternal phase drifts between gates [ljj, |2fJ, |21| . 

At the heart of Grover's algorithm is the "oracle 
query," which quickly checks if a proposed input "x" 
is a solution to the search problem. The oracle marks 
a particular component of a quantum superposition by 
flipping the sign of its amplitude. Following the ora- 
cle, a number of quantum operations amplify the weight- 
ing of the marked state independent of which state is 



FIG. 1: Schematic diagram of Grover's quantum search al- 
gorithm over a space of n qubits (JV = 2™ entries). Initially, 
all qubits are prepared in the state |000...0). (a) A global 
Hadamard gate prepares an equal superposition of all states, 
(b) The oracle (shaded in yellow) flips the sign of the am- 
plitude corresponding to the marked element, represented by 
the n— bit binary number x. (c) Two global Hadamard gates 
surround an additional phase gate (shaded in red) that flips 
the sign of the amplitude corresponding to the initial state 
|000...0), amplifying the weight of the marked state. Steps 
(b) and (c) are repeated in sequence a prescribed number of 
times, and finally (d) the qubits are measured. An example 
of the distributions of quantum amplitudes at each stage are 
depicted at the right. 



marked (see Fig. After many iterations of this 

query/amplification process, the marked state accumu- 
lates nearly all of the weight and is revealed following 
a measurement. The required number of queries is the 
integer closest to 4sm _ 1 ^ Af _ 1/2) -^ [H For JV > 1, the 
marked element would thus appear with high probability 
after approximately ny/N /4 iterations, and for the spe- 
cial case of N = 4 elements, a single query would provide 
the marked element with unit probability. Classically, a 
single query of a 4-element search space followed by a 
guess can only result in a successful outcome with 50% 
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probability. 

A standard quantum circuit for the Grover search al- 
gorithm for N=4 entries is shown in FigureEJa) ,1]. This 
scheme uses a third ancilla bit which marks one of the 
database elements through a Toffoli gate that effectively 
flips the sign of the marked element if and only if the two 
bit input is a solution to the problem. The oracle scheme 
to mark each of the four possibilities is shown below the 
circuit. The remainder of the circuit (shaded in dark red) 
amplifies the weighting of the marked state, with the op- 
erations between the Hadamard gates flipping the sign of 
the amplitude of the 1 00) state. 

Fig. EJR) shows the experimental implementation of 
the algorithm for N=4 search elements. The M0lmer- 
S0rensen (M-S) entangling gate [23| is adapted to the 
algorithm for both the marking and state amplification 
steps and is accompanied by a variety of single qubit ro- 
tations. This circuit is identical in function to the circuit 
of Fig. |2t a) but does not use an ancilla qubit [24j, First, 
each qubit is initialized to the |0) state through optical 
pumping techniques. Next, an equal superposition of all 
qubit states is prepared with a global single qubit ro- 
tation of 7r/2, replacing the Hadamard gates of Figs 
and |3[a) . The oracle function (shaded in yellow) then 
marks one of the four possible states \xx) by flipping its 
sign. This is accomplished with a controlled-z phase gate 
(shaded in blue) containing the M-S gate, that flips the 
sign of the |11) state. The phase gate is surrounded by 
rotations (shaded in gray) that swap the marked state 
\xx) with the |11) state, depending on the settings of the 
of the rotation angles depicted in the table at the bottom 
of Fig l^b) . After the oracle query, a global rotation of 
7r/2 (with a phase of -it/ 2 relative to the initial rotation) 
followed by a second application of the M-S gate ampli- 
fies the weighting of the marked state \xx) (operations 
shaded in red). 

The experiment is performed with two 111 Cd + ions 
confined in a three-layer linear ion trap with axial fre- 
quency w z /27r=2.0MHz 0110. The S 1/2 ground 
state hyperfine levels \F — 0,m/ = 0) (denoted by |0)) 
and \F = l,m/ = 0) (denoted by |1)), separated in fre- 
quency by ujo/2tt= 14.5 GHz, serve as qubit levels. These 
"clock" states are relatively insensitive to magnetic field 
noise, obviating the need for "spin echo" pulses in the 
implementation of the algorithm [l^. A variety of op- 
tical and microwave sources are used to control, entan- 
gle, and measure the qubits. In particular, pairs of off- 
resonant laser beams with a beatnote near luq drive stim- 
ulated Raman transitions that couple the qubits to their 
collective motion and give rise to entangling quantum 
gates jl3,|2^|. Both the center-of-mass and stretch modes 
of axial motion are initialized to near the ground state 
via stimulated Raman cooling, with average vibrational 
numbers of n=0.28 and 0.12, respectively j^. Single 
qubit rotations are accomplished through resonant mi- 
crowave pulses (tuned to near loq) and composite pulse 
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FIG. 2: Quantum circuit to implement Grover's searching 
algorithm for N=4 entries [l]. (a.) Theoretical circuit us- 
ing a third ancilla bit and standard gates including the Had- 
damard gate (H), the generalized Toffoli gate, a bit flip,X, 
and a controlled-NOT gate. The Toffoli gate implements the 
oracle (shaded in yellow), where the scheme to mark each of 
the four possibilities is shown below the circuit. The remain- 
der of the circuit (shaded in red) amplifies the weighting of 
the marked state, (b.) The experimental circuit to implement 
the algorithm for n=2 qubits, where R(#, (f>) is a rotation on 
the Bloch sphere |32||. R z {4>) is a phase rotation about the z- 
axis, and Gms is the M0lmer-S0rensen entangling gate. The 
yellow shaded box identifies the oracle, where the value of the 
variables a and j3 (given in the table), determine which state 
is marked |33|. The remainder of the circuit (shaded in red) 
amplifies the weighting of the marked state. 



sequences involving a single tightly-focused off-resonant 
laser. The microwaves are applied to both qubits simul- 
taneously and provide global single qubit rotations. The 
off-resonant laser, which addresses one ion more strongly 
than the other, creates a differential AC Stark shift be- 
tween the two qubits, imparting a relative phase shift of 
7T between the two qubits. Combined with appropriate 
microwave pulses, this allows arbitrary single qubit rota- 
tions. Finally, each qubit is detected with greater than 
97% efficiency by uniformly illuminating the ions with 
resonant laser radiation and observing the ion fluores- 
cence on an intensified CCD camera. 

The M0lmer-S0rensen gate directly entangles the clock 
state qubits and is insensitive to the relative optical phase 
of the Raman laser beams between gates 19, 20, 26] . This 
is an important consideration when multiple entangling 
gates are implemented because it suppresses decoherence 
from magnetic fields and optical phase noise that may 
fluctuate from gate to gate. The M-S entangling gate is 
realized by applying multiple sets of Raman beatnotes to 
the ions, simultaneously driving the first lower and upper 
motional sidebands for a particular duration. This entan- 
gles the spin states via their collective motional mode, in 
this case the stretch mode [23| . The evolution of the four 
basis states is given by: 
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The desired states are produced with a fidelity of approx- 
imately 80%. There are additional phases not included in 
the above equation pertaining to the ion-ion spacing, the 
phase of the oscillator that defines the Raman beam beat- 
note, and Stark shifts from the applied Raman beams 
[l9l l21j . These phases are set to zero for the present case 
since they are fixed before the experiment is run by first 
synchronizing the phase of the entangling gate with the 
phase of microwave 7r/2 pulses [l^. Then the phases of 
the two entangling gates used in the experiment are syn- 
chronized to each other through a Ramsey experiment. 
Note that phase noise from the Raman beam path fluc- 
tuations is supressed due to the spectral arrangement of 
the Raman sidebands [i3,[2^|. The timescale for each op- 
eration in the algorithm is as follows: 10/is for a global 
microwave rotation, 20^ts for a differential single qubit 
rotation, and 140/xs for the M0lmer-S0rensen two qubit 
entangling gate, giving a total of ~ 380^s to complete 
the 20 pulses that form the algorithm. 

There are several approaches to gauging the perfor- 
mance of the algorithm implementation. One method 
is to compare the algorithm's success at recovering the 
marked state with the best that can be achieved classi- 
cally. The classical counterpart is a simple shell game: 
suppose a marble is hidden under one of four shells, and 
after a single query the location of the marble is guessed. 
Under these conditions, the best classical approach gives 
an average probability of success P c i = 1/4 + 3/4(1/3) = 
0.50, because 1/4 of the time the query will give the cor- 
rect location of the marble while 3/4 of the time a guess 
must be made amongst the three remaining choices each 
with 1/3 probability of choosing the correct location. If 
Grover's algorithm is used, the answer to the single query 
would result in a 100% success rate at 'guessing' the mar- 
ble's location. As can be seen in figure EJa) the marked 
state is recovered with an averaged probability over the 
four markings of 60(2)%, surpassing the classical limit of 
50%. 

It is interesting to consider the output of the algorithm 
when the final entangling gate used for state amplifica- 
tion is omitted. This situation shows how well the al- 
gorithm can do with only single qubit rotations outside 
the oracle. This scenario lies between the classical and 
quantum searches described above since entanglement is 
not used outside the oracle but quantum superpositions 
are used to find the marked element. In this case it can 
be shown that quantum mechanics without entanglement 
can do no better than what can be achieved with clas- 
sical means: both methods have the outcome of finding 
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FIG. 3: (a.) Output of the algorithm. The conditional 
probability of measuring each of four output states given 
one was marked is shown in sequence from top to bottom 
|00), |01), 1 10), |11). Each of the four data sets shows the 
distribution of measurements averaged over 500 trials. The 
marked state was recovered on average with 60(2)%, com- 
pared to unit probability for the ideal quantum algorithm and 
50% for the best possible classical algorithm, (b.) Output of 
the algorithm without the final entangling gate. This shows 
that the fidelity of the oracle is about 80%. Each of the four 
data sets was also averaged over 500 trials. The experimental 
average to recover the marked state is 41(2)% with the theo- 
retical limit of 50%, both of which are less than the 60% from 
(a). The quoted errors are statistical. 



the marked state with only 50% probability, assuming 
a perfect oracle. In addition, this diagnostic allows the 
performance of the oracle itself to be characterized. The 
rotations following the oracle convert the marked state 
into one of four Bell states each of which yields a max- 
imum probability of 50% to recover the marked state. 
Figure |3b shows that the marked state is recovered with 
an average of 42(1)% probability, implying the oracle it- 
self has a fidelity of roughly 80%. 

The above figures of merit focus on the mean success 
probability and neglect the information content inherent 
in the distributions of figure^. The mutual information 
between the marking of the state and the measurement 
can be used to characterize this correlation and hence is 
another measure of the a lgor ithm's success p| (for other 
figures of merit see ref. |24|). The mutual information 
measures how much information two random variables, x, 
the measurement, and y, the marking, have in common. 
It is defined by: H(x:y)=ff (cc) + H{y) — H(x,y), where 
H(x, y)=-J2 x ,y P( x >y) lo S2 P( x >y) is the joint Shannon en- 
tropy between the two distributions, p(x, y) — p{x)p{y\x) 
is the joint probability distribution of x and y, and p{y\x) 
is the conditional probability of y having been marked 
given that x was measured. H (^—-J^xPi^) 1°§2 p(x) 
and H(y)=-^2 y p(y) log2 p(y) are the Shannon entropies 
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of the individual variables. Classically the mutual in- 
formation acquired after a single query of the oracle 
is H(x:y)=0.251og 2 (0.25)-0.751og 2 (0.75)=0.81 bits, mean- 
ing, on average, 0.81 bits of information are gained upon 
measurement. The ideal quantum algorithm would yield 
two bits of information upon measurement. For the data 
in figOUi the mutual information is 0.44, so on average 
only about a half a bit of information is gained. Even 
though less information is gained per measurement than 
the classical case, the probability of finding the marked 
state in the experiment still exceeds the classical limit. 

As the data in fig. [3^ shows, the marked state is not 
recovered with unit probability. A large part of this in- 
fidelity is due to the M0lmer-S0rensen entangling gate. 
Each instance of the M0lmer-S0rensen gate has a fidelity 
of about 80%, and since there are two such gates in the 
algorithm, overall fidelities of approximately 60% are ex- 
pected. The main sources of decoherence during the gate 
are spontaneous emission from off-resonant coupling to 
the excited state and fluctuating AC Stark shifts from the 
Raman beams that drive the entangling gateji^. Both of 
these induced decoherence sources can be suppressed by 
increasing the detuning of the Raman beams from the 
excited state, at the expense of slowing the gate. We 
choose the detuning to strike a balance between these in- 
duced decoherence sources and other slowly varying noise 
sources, such as motional heating [2^, [23, fluctuating 
magnetic fields, and microwave oscillator phase drifts. 
Additional power in the Raman laser beams accompa- 
nied by larger detunings could suppress decoherence from 
spontaneous emission and AC Stark shifts while main- 
taining a reasonable gate speed (see [20J for more details). 
Fluctuating AC Stark shifts during the differential single 
qubit rotations due to technical intensity fluctuations and 
beam pointing instabilities add infidelities to the experi- 
ment on the order of 5-10%. 

For Grover's algorithm to be useful it needs to extend 
beyond a few qubits. Using a quantum circuit similar to 
figGK- an fi-qubit Grover algorithm can be implemented 
with n-qubit Toffoli gates, a series of two qubit gates, and 
single qubit rotations. It has been shown that an n-qubit 
Toffoli gate can be constructed with single qubit gates 
and controlled-NOT gates with order n basic operations 
[2^| . A controlled-NOT gate can be constructed from the 
M-S entangling gate through the following sequence [2^ : 
[i? 2 (7r/2, 0), i?i(7r/2, tt), R 2 (tt/2, tt), G M s, Ri(n/2, 0),* 
R 2 (n/2, 0), R z1 {-k/2),R z2 {-k/2), R 2 (n/2, -ir),R zl {n)}, 
where ^4=1,2(0, (f>) is a rotation of ion i by angle 9 and 
phase <j), R z i (4>) is a z- rotation of ion i by angle </>, and 
Gms is the M0lmer-S0rensen entangling gate. Since 
the ion system is scalable to a large number of qubits 
it is feasible to construct an efficient n-qubit Grover 
algorithm where each iteration scales polynomially with 
n. In this case, the isolation of individual ions could be 
accomplished through tight focusing of laser beams j2j| 



or the shuttling of ions between separated trap zones 
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Under this convention the basis states evolve as: 
|0> -> cos(0/2)|O) + e^sm(0/2)|l) and |1) 
-e-^sin{6 /2)]0) + cos(0/2)|l) 

In the experiment the first set of rotations (shaded in 
gray) are accompanied by extra z-rotations. These addi- 
tional phases do not affect the outcome of the oracle and 
are suppressed in the figure for simplicity. 



